Dynamic simulation of hydrodynamically interacting particles

A general method for computing the hydrodynamic interactions among N suspended particles, under the condition of vanishingly small particle Reynolds number, is presented. The method accounts for both near-field lubrication effects and the dominant many-body interactions. The many-body hydrodynamic interactions reproduce the screening characteristic of porous media and the ‘effective viscosity ’ of free suspensions. The method is accurate and computationally efficient, permitting the dynamic simulation of arbitrarily configured many-particle systems. The hydrodynamic interactions calculated are shown to agree well with available exact calculations for small numbers of particles and to reproduce slender-body theory for linear chains of particles. The method can be used to determine static (i.e. configuration specific) and dynamic properties of suspended particles that interact through both hydrodynamic and non-hydrodynamic forces, where the latter may be any type of Brownian, colloidal, interparticle or external force. The method is also readily extended to dynamically simulate both unbounded and bounded suspensions.